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We develop and implement a model for decoherence in time-dependent transport. Inspired in a 
dynamical formulation of the Landauer-Biittiker equations, it boils down into a form of wave function 
that undergoes a smooth stochastic drift of the phase in a local basis, the quantum-drift (QD) model. 

This drift is nothing else but a local energy fluctuation. Unlike quantum-jumps (QJ) models, no 
jumps are present in the density as the evolution is unitary. As a first application, we address the 
transport through a resonant state 10) that undergoes decoherence. Its numerical resolution shows 
the equivalence with the decoherent steady-state transport in presence of a Biittiker’s voltage probe. 

In order to test the dynamics we consider two many-spin systems, which are cases of experimental 
interest, where a local energy fluctuation is a natural phenomenon. A two-spin system is reduced to a 
two-level system (TLS) that oscillates among 10) = |j\|_) and |1) = |JJ). We show that the QD model 
recovers not only the exponential damping of the oscillations in the low perturbation regime, but 
also the non-trivial bifurcation of the damping rates at a critical point, i.e., the quantum dynamical 
phase transition. We also address the spin-wave-like dynamics of local polarization in a spin chain. 

By averaging over N s realizations, the QD solution has about half the dispersion respect to the 
mean dynamics than QJ. By evaluating the Loschmidt echo (LE), we find that the pure states |0) 

= |j\f) and 11) = |4A) are quite robust against the local decoherence. In contrast, the LE, and 
hence coherence, decays faster when the system is in a superposition state (|j\f) i Iff)) /\/2, which 
is consistent with the general trend recently observed in spin systems through NMR. Because of its 
simple implementation, the method is well suited to assess decoherent transport problems as well 
as to include decoherence in both one-body and many-body dynamics. 

PACS numbers: 03.65.Yz, 72.10.—d, 42.50.Lc 


I. INTRODUCTION 

The last decade has seen an increasing demand to de¬ 
scribe quantum dynamics on a variety of complex sys¬ 
tems in the presence of an environment. Among them 
are atomic systems in optical lattices EES , networks of 
interacting spins [3] and charge and magnetization dy¬ 
namics of nanoscopic devices in a transport set up 0-0 ■ 

The most common way to deal with environmental de¬ 
coherence in small closed systems, is the master equation 
for the density matrix in a Lindblad form (3, 0 • In or¬ 
der to deal with bigger systems, the most standard ap¬ 
proaches implement the Redfield theory for a relaxation 
superoperator which does not ensure the strict unitarity 
of the Lindbad form 00. Alternatively, some works 
pointed to strategies for computing the evolution of an 
open system based on the stochastic dynamics of a state 
vector that suffers instantaneous quantum jumps (QJ) 
[M3- Indeed, for large systems, it was shown that 
the stochastic method is faster than the density matrix 
implementations 0 - 

With regard to decoherent steady state transport, the 
traditional evaluation in terms of the Kubo linear re¬ 
sponse 00 was progressively replaced by the Lan- 
dauer’s motto that conductance is a quantum transmit- 
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tance fl9 ] . This is implemented in the Landauer-Biittiker 
(LB) scattering formulation where decoherent processes 
are induced by current conserving voltmeters. These im¬ 
pose a self-consistent reinjection of the electrons that 
ensures the current cancellation. Biittiker had the idea 
that this reinjection describes a decoherent process 0 . 
a concept extended by D’Amato and Pastawski to deal 
with more general situations 0,[0 becoming a popular 
tool to address decoherent transport 0-0 ■ A dynami¬ 
cal formulation, called the generalized Landauer-Biittiker 
equation (GLBE), based on the Kadanoff-Baym-Keldysh 
(KBK) quantum field theory for non-equilibrium pro¬ 
cesses was then developed [0. The GLBE seeks to find, 
in a linear response, the non-equilibrium Keldysh density 
function which is proportional to the density matrix. In 
recent years, there has been a burst of progress in the 
use of the KBK formulation of non-equilibrium dynam¬ 
ics in a framework consistent with ab initio calculations 
0,0. Nonetheless, numerical solutions of many-body 
systems become excessively demanding. In particular, 
they involve time integrals of self-consistent memory ker¬ 
nels 0-0| . Additionally, in strongly interacting many- 
body systems, which are beyond a mean-field description, 
such as spin systems, it would involve costly averages over 
the participating configurations. 

The above limitations can be overcome by resorting 
to two strategies. The Trotter-Suzuki step-by-step evo¬ 
lution makes practical some calculations of the Keldysh 
density function 0, 0. The ensemble average is by- 
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passed by a recent strategy, dubbed “quantum paral¬ 
lelism” , which uses a single wave function that is a super¬ 
position of all the states participating in the statistical 
ensemble a However, these strategies are limited to 
coherent dynamics. 

In this work, we propose a stochastic model that 
extends the Biittiker-D’Amato-Pastawski approach to 
evaluate decoherent time-dependent problems. Our 
quantum-drift (QD) model is based on a wave func¬ 
tion stochastic dynamics. Within a discrete-time set up, 
we impose an incoherent reinjection that ensures den¬ 
sity conservation. Then, as in quantum parallelism, we 
propose that the wave function should sum up both a 
coherent and an incoherent part. Thus, being fairly rep¬ 
resentative of a set of stochastic interaction “histories”, 
this wave function does not present jumps but smooth 
drifts in a single unitary evolution. Further averages over 
realizations are performed only in the needed amount ac¬ 
cording to the addressed observable. 

This paper is organized as follows. In Sec. II we re¬ 
view the basis of the Biittiker’s model for transport in 
open systems. In Sec. Ill we present the basis of our QD 
model. By addressing a wave packet dynamics through 
a double barrier resonant tunneling device (DBRTD), in 
Sec. IV, we show that QD recovers the decoherent steady 
state transmittances of the Biittiker formulation. In Sec. 

V we compare the QD dynamics with the Keldysh so¬ 
lution in a decoherent closed system with a simple but 
non-trivial situation: a two-level system (TLS) that un¬ 
dergoes a quantum dynamical phase transition (QDPT) 
[3^ |. In Sec. VI, we perform a many-spin calculation for 
a case of experimental and theoretical interest, that of a 
quantum channel, showing the agreement of QD and QJs 
and comparing their numerical performance. Finally, in 
Sec. VII we use the QD model to evaluate decoherence 
through the Loschmidt echo (LE) in the TLS of the Sec. 

V under local decoherence processes 0. This allows us 
to show that while Rabi oscillations decay uniformly, de¬ 
coherence is not uniform: it affects the system only while 
it goes through non-local superposition. We present the 
final discussion in Sec. VII. 


II. DECOHERENT QUANTUM TRANSPORT: 

REINJECTION, PARALLELISM, 
ATTENUATION, AND ENERGY BROADENING. 

The first phenomenological model for decoherence was 
developed for Biittiker in the context of electronic trans¬ 
port in phase-coherent mesoscopic systems fl9j . He re¬ 
alized that decoherence may be introduced by including 
a terminal connected to a voltmeter. This key idea may 
be readily visualized by considering a three-terminal cir¬ 
cuit. There, two terminals (source L and drain R) pro¬ 
vide an infinitesimal voltage difference that produce a 
current through the system while the third one, is con¬ 
nected to the voltmeter. The electrons from the left ( L ) 
and right ( R ) leads that enter into the voltmeter undergo 


a decoherent process. An appropriate local chemical po¬ 
tential at <p (i.e. the measured voltage) ensures current 
cancellation, i.e. a reinjected electron compensates each 
electron that flies into the voltmeter. Then, this electron 
does not keep any memory or phase correlation with re¬ 
spect to the electrons inside the sample. 

Let consider that represents the quantum coherent 
transmittance from the j to i channel connecting them to 
the reservoirs ( i 7 ^ j, take the values L , R or </>). The ap¬ 
plication of the Landauer-Biittiker equations for a system 
with one voltmeter results in a transmittance through the 
system given by 


Trl 


Tr^T^l 


coherent 

decoherent 


Tr,/, + T^l 


Trl + 


(1) 


The first term is the probability that the particle travels 
from L to R without undergoing a decoherent process 
at the voltmeter </>. The second term accounts for those 
electrons that have interacted with the environment at 
cj). Indeed, one can recognize it as the conductance of 
two parallel pathways. One of them with a conductance 
( 2 e/h)TRL, while the other one adds a series of two con¬ 
ductances, (2e/h)TR C j > and ( 2e/h)T ( f > L■ 


These results adopt a more concrete form by introduc¬ 
ing an explicit model. Let us consider a quantum dot 
where E 0 is the relevant eigenstate energy, e the Fermi 
energy, Eq the local energy properly shifted by the pres¬ 
ence of the contacts, and and T#, the energy un¬ 
certainties produced by the escape toward the left and 
right leads, respectively. The dot’s Green’s function is de¬ 
fined in terms of the effective Hamiltonian parameters as: 

G(e) = [(e — Eq) + i(To + T^)] , where To = 

is the natural width due to the presence of the leads 
[H, 0. Thus, each of the transmittances used above 
may be written explicitly in terms of G{e) by using the 
Fisher-Lee formula Tij{e) = 2Fj \G(e)\ 2 2Tj with i 7 ^ j. 


If Tftl is the transmittance from L to R that accounts for 
electrons in absence of the voltmeter (i.e. when T^ = 0 ), 
the coherent part Trl in Eq. [T]can be further written as 
the product of and an attenuation factor (1 — A(e)) 
f3E| . Thus, the effective transmittance is written in terms 
of a attenuated coherent part plus an incoherent one. 

D’Amato and Pastawski’s (DP) model generalize these 
ideas and introduce an effective Hamiltonian that consti¬ 
tutes a microscopic model for decoherence in the steady 
state [II|. In this case, the isolated system is described 
by a Hamiltonian Ho. DP identify the escape of the elec¬ 
trons towards the voltmeter with their interaction with 
the infinite degrees of freedom of an environment. Here, 
decoherence is induced by local processes (e.g., a volt¬ 
age probe, a local phonon bath) in the Fermi golden rule 
(FGR) approximation. These interactions produce an 
energy uncertainty F^ for each local state with a rate of 
system-environment interaction, 1/tre = ‘2T ( j > /h : which 
has an irreversible character |2l|, [26, 341. 
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The local density of states (LDoS) is calculated from 
the dot’s Green’s function 


iV 0 (e) 


- 

7 r 


1 _To 

n(e-E 0 )i + r 2 0 


( 2 ) 


By including the system-environment interaction, the 
LDoS acquire an extra energy uncertainty or broadening 
Then, the LDoS in presence of decoherence, No{e), 
is obtained from -ZVo(e) by replacing the characteristic 
width, To, byTo+T^. On the other hand, Nq may be 
obtained by considering that individual decoherent pro¬ 
cesses shifts the resonances in an amount A E from Eq. 
Then, by averaging over the possible A E, 

N 0 (e) =Jn 0 {£- AE)P(AE)dAE, (3) 
with P(AE) = - -^-. (4) 

tt(A £) 2 + V’ 


where, P{AE) is a Lorentzian probability distribution 
for the shifts A E. Similar broadening occurs in other 
observables that depend on energy, such us correlation 
functions. 


III. THE QUANTUM DRIFT MODEL 

By using the Trotter-Suzuki expansion, the quantum 
dynamics is obtained by the sequential application of uni¬ 
tary evolution operators that transform the initial state 
in small time-steps, dt. If the system-environment in¬ 
teraction has a rate 1/tse, during each interval dt, the 
particle has a probability p = d t/rsE to undergo a de¬ 
coherent process and a probability (1 — p) to survive it 
mm. Let us consider a single state |0) that may un¬ 
dergo a decoherent process. Thus, after a time dt, the 
coherent amplitude is reduced by the factor a/1 — p owing 
to decoherent processes. In order to conserve the density, 
and in accordance with the Landauer-Biittiker picture, 
the wave function must include a term that accounts for 
the decoherent reinjection. Thus it has a random phase 0 
drawn from some distribution Pg. We can represent both 
the coherent and the incoherent contributions in the same 
wave function, resembling quantum parallelism. This is, 


for some random phase /3o- 

Notice that, in a Trotter-Suzuki evolution, the phase 
shift e l/3 ° is actually a correction AEq = hfio/dt in the 
energy of the state |0). Equation [3] shows that the sin¬ 
gle level coupled to an environment acquires an energy 
uncertainty T^, which in a FGR approximation, is char¬ 
acterized by a Lorentzian shape. This, in turn, is as¬ 
similated to a distribution of instantaneous energy shifts 
AEq drawn from the Lorentzian distribution. In our 
model, the correction AEq is taken to be a random num¬ 
ber that varies step-by-step to represent the uncertainty 
introduced by the environment. Thus, the probability 
distribution Pp 0 is 


pm = - 

7r 


r^d t/h 

$ + (r*dW 


( 8 ) 


Therefore, the key decoherent processes are the highly 
improbable processes that involve a large A Eo (the tails 
of the Lorentzian). 

This proposal can be extended to all the levels E n 
of a Hamiltonian in an arbitrary basis. In particu¬ 
lar, in a tight-binding basis, each site energy E n ac¬ 
quires a Lorentzian energy uncertainty T^ n and these 
are perturbed with a random energy A E n . More for¬ 
mally, we can define E as a diagonal operator where 
2n,n' = AE n 5 n ^ n i. For a N x N matrix Hamiltonian 
Hq we consider an effective instantaneous Hamiltonian 
H e ff . = Hq + E. Thus, we obtain the unitary evolution 
operator in a Trotter-Suzuki expansion, 

U(dt) = e - ii? e//dtA ( 9 ) 

~ e -tidt/h e -iH 0 dt/h = /7 s (di)U 0 (df). (10) 

We can define the decoherence operator as Z7 e = 
exp[—iEdf//i], that is unitary, it conserves the density 
and, thus, the density does not presents jumps. However, 
observables involving two-sites correlation (e.g. currents 
and momentum) do have jumps. However, these last are 
smoothed out by taking the ensemble average. 

In summary, the prescription to include decoherence 
in a quantum dynamics is to include in every time step a 
random correction, /?„, to the phase of each local state. 
This correction has a distribution probability given by 
the P/ 3 n of Eq. [5] Thus, the evolution of a wave function 
is performed by 


— coh. incoh. 

W = W + V>o (5) 

= (a / 1 -P + x e e ' e ) i>o ( 6 ) 

where ipo = (0| The cross terms cancellation in the 

ensemble average is ensured by f Pg^fl — p\ge l6 dO = 0. 
In any case, the coefficient A g should be chosen to account 
for density conservation |-/l — p + \ge l9 \ = 1. Thus, 

V’o = e 10o tp o, (7) 


N t 


= II' 


,—iTidt/h—iHodt/h 


WO)}, (if) 


j=i 


where N t = t/dt. 


IV. DECOHERENT TRANSPORT: 
D’AMATO-PASTAWSKI TRANSMITTANCE 

Let us first test our model in the system that in¬ 
spired it: decoherent transport through a double barrier 
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resonant tunneling device (DBRTD). In a tight-binding 
scheme, this is represented by one resonant site of energy 
Eq = 0 coupled to two semi-infinite leads of bandwidth 
4V (where V is the unit of energy), L and R that act 
as current source and drain. The tunneling amplitudes 
through the barriers are Vl and Vr. Thus, the tight- 
binding Hamiltonian is 


Hq — EqCq Cq V^c+iCq + c_i) ^(c+co + cjcj (12) 


OO 

- J2 v (c+ +1 c n + c+c n+1 ) - 

n =1 


— OO 

E 

n=— 1 


V(£+C n + C+C n+1 ). 



e (units of V) 


To evaluate the transmittance from a dynamical cal¬ 
culation, we build a Gaussian wave packet well inside 
the left lead. A wide wave packet ensures a well defined 
energy. The transmission coefficient is obtained by inte¬ 
grating the density at the right side after the wave packet 
has been transmitted or reflected. This transmittance is 
equivalent to the steady-state analytic result of the Fisher 
and Lee formula [34}. 

Decoherence is introduced only at the resonant level as 
described in the previous section during the whole evo¬ 
lution. In Fig. |T]we compare the QD results with those 
resulting from the Biittiker’s solution of Eq. |T] We plot 
these quantities for different decoherence strengths T^. 
These are in a very good agreement, made even more 
valuable by considering that the number of realizations 
in the average was of the order of 10. This is because the 
observable of transmitted density involves a spatial inte¬ 
gration. The same self-average is observed for the deco- 
herent conductance in long one-dimensional (ID) wires. 
In that case even a single realization is enough to repro¬ 
duce the known results in this problem pH ]. 

The QD method fits the theoretical values to the de¬ 
sired precision. The only difference that one might no¬ 
tice in certain specific cases, such as narrow peaks, would 
arise from the fact that wave packets are built from states 
within an energy range as a consequence of the uncer¬ 
tainty principle. On the other hand, scattering theory, 
by using asymptotic plane waves, has no energy uncer¬ 
tainty. This example is representative of a wide variety 
of steady state problems that can be solved with the QD 
method, at only the cost of numerical resources. Indeed, 
we tested QD on extended systems, a situation where an 
implementation of the QJ would impose excessive local 
fluctuations. We found that in QD, much as in quan¬ 
tum parallelism [29}, the collective observables have a 
tendency to self-average. This makes our model a very 
promising tool to evaluate decoherent dynamics in ex¬ 
tended systems and in many-body problems. However, 
the true advantage of QD starts to be appreciated in ad¬ 
dressing time-dependent problems, as we do in the next 
section. 


FIG. 1. (Color online) Decoherent transmittances obtained 
for two values of the decoherence strength. For T^/V = 0.01, 
the Biittiker transmittance is plotted with the short-dot line 
and the transmittance of the QD model is plotted with the red 
triangles for a number Ns = 5 of realizations in the average. 
For T^/V = 0.3, the Biittiker transmittance is showed by 
the solid line and the QD result, with the blue circles for a 
number of Ns = 50 realizations. The other parameters are: 
V = 1, E 0 = 0, V L = Vr = 0.15, E r = 0 


V. QUANTUM DYNAMICAL PHASE 
TRANSITION IN A TWO LEVEL SYSTEM. 

Let us consider a two-level system (TLS) that de¬ 
scribes charge or spin dynamics, [1, 30] say states |0) 
— It'i') = c+ |' vacuum) and |1) = |4-t) = \vacuum), 

with degenerate energy Eq and an interaction V mixing 
them. Such simple system was seen to have a non-trivial 
dynamics when one of its levels interacts with an envi¬ 
ronment of spins: a quantum dynamical phase transition 
(QDPT) [32]. In a QDPT certain observables present a 
non-analytic dependence on the system-environment in¬ 
teraction strength. The QDPT was missed in a solution 
for the density matrix in the usual secular approxima¬ 
tion of the Redheld theory [3(| but showed up in a QJ 
variant [32]. When the TLS suffers the asymmetric in¬ 
teraction of an environment, the QDPT already appears 
in the spectrum of the effective non-Hermitian Hamilto¬ 
nian [371 ]. However, if the interaction is symmetrical, the 
QDPT only occurs in the density matrix if positivity is 
ensured [38}. Thus, obtaining the QDPT in a model with 
symmetric interaction with the environment constitutes 
a definitive test for the QD method. 

The Hamiltonian of the TLS is 

H 0 = Eq(cq c 0 + cfcj - U(c^c 0 + c^cj. (13) 

The survival probability of an excitation with an initial 
state \ip (0)) = |0), i.e. the diagonal element of the den- 
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sity matrix, is 


Poo (t) 


( 0 | e -' 1%t/h | 0 ) , 


= 1/2 + 1/2 coscoot. 


(14) 

(15) 


Here, we can observe that Poo(t) has Rabi oscillations 
HH with frequency co o = 2 V/H and period T = irh/V. 

Let us consider now an environment that interacts in¬ 
dependently with each state with a rate described by the 
Fermi golden rule 1 /tse = 2r ^/H, where, F^ is the en¬ 
ergy uncertainty of the level. Physical implications of 
this model are discussed in the next section in the con¬ 
text of a spin system. The numerical evolution of the 
TLS is now performed by choosing as the width of 
the Lorentzian distribution. 

Here we will compare our QD method, where the de- 
coherent survival probability is 


Poo (0 


N t 

(0| e-ifi-dt/ftg-iSodt/Rio) 

n =1 


2 


(16) 


with the analytic solution of the GLBE. This last was 
analytically solved for this problem in Refs. [22j. ItE [ii 8 j . 
giving for the survival probability 


Poo(t) 


A _l ip-rv/ft 
2 2 


cos (cot) -— sin (cot) 

2 co 


(17) 


Thus, the oscillations of both the diagonal and non¬ 
diagonal elements of the density matrix oscillate with a 
frequency co, which is lower than the Rabi frequency coo 
according to 


co = \Jcol- (T^/nf . (18) 

This evidences that the oscillation frequency of a TLS 
exhibits a non-analytic behavior. The frequency co takes 
real values provided that T^/H < co 0 (underdamped 
regime). Beyond this value, i.e., for T^/H > co 0 (over¬ 
damped regime), Re(co) = 0 and thus the oscillations are 
fully suppressed, and Poo(t) is the sum of two exponential 
decays: 


*"(‘>- 5 - 5 ^ 


-lit 


7i 


0-121 


2(71 - 72 ) 

where the decay rate 71 ( 2 ) is 

71 ( 2 ) = jc (t* ± - (hoo) 2 ) ■ 


, (19) 


( 20 ) 


Note that, at short times, Poo(t) is always of the form 
1 — co 2 t 2 /4 = 1 — V 2 t 2 /h 2 , which is characteristic of a 
quantum evolution without perturbations. This is be¬ 
cause, at short times, the environment interplay has a 
small cumulative effect on the survival probability, which 


is still determined by the unperturbed quantum dynam¬ 
ics. In a strongly decoherent regime, T^/h co 0 , the 
decay rates tend to 71 ~ 2T^/h, 72 — hco 2 /2T ( j 1 , defining 
a short-time decay rate, 71 , and a rate 72 , that dominates 
the long times as Poo(t) oc e _72t . Both exponential terms 
are needed to obtain the whole evolution. An equivalent 
solution for the QDPT may be obtained by considering 
models for environmental noise j4o| . 



FIG. 2. (Color online) (a) Theoretical rates of decay pre¬ 
dicted by the Eqs. 1171 and l20l Red points indicate the rates 
of the GLBE solution for (b) — (d). The QDPT occurs as 
a bifurcation in the rates for the critical value r^/F = 2. 
The dotted line represents the asymptotic value 2T ,/,/h. In 
( b) — (d), we show the survival probabilities Poo(t) (thick blue 
line), the GLBE solutions (thin black line), and individual 
realizations (thin gray line) in the underdamped regime, in 
( b ) with F^/F = 0.1 and (c) with F^/F = 1, and in the over¬ 
damped regime, (d) with T^/V = 2.1. Individual realizations 
tend to preserve the oscillations. These do not have jumps 
in the densities but in the slopes. The ensemble average is 
taken over Ns = 100 realizations. For Poo(t) ~ 1/2 are clear 
the typical fluctuations of the order of 1 /y/Ns- The average 
Poo(t) tends to the GLBE solution by increasing Ns- 

In Fig. m.a) the decay rates are shown. The QDPT is 
manifested as the bifurcation in these rates. The fitted 
rates are shown as points superposed to the theoretical 
curve. In Fig. [2j £»), we show the Rabi oscillations of the 
average survival probability, PooWi which are exponen¬ 
tially attenuated with r _1 ~ T^/h = 1/2 (tse) -1 - This 
is the most common example of decoherence in TLS’s. 
The Poo(t) and the fitted decay rates match perfectly the 
GLBE solution. A single realization of QD method is also 
shown. Notice that there are no “jumps” in the survival 
probability, and an oscillatory-like behavior dominates 
the whole evolution. In Fig. [3 c) we show the Poo(t) i n 
the underdamped regime for a value of = IV, where, 
the oscillatory behavior is small. In Fig. (3^) we show 
Poo(t) in the overdamped regime, for = 2.IF. We 
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can identify the initial quadratic behavior and, at large 
times, the exponential decay with the rate 72 oc l/T^. 
The larger T^, the slower Poo(t) decays. This is a sig¬ 
nature of the quantum Zeno effect in which the system 
is continuously perturbed freezing its evolution close to 
the initial condition. By increasing T^, each single real¬ 
ization seems to be a stochastic process while preserving 
the quadratic initial starting. Single realizations do not 
present jumps in the density but in the correlations, seen 
as sudden changes on the slopes. 

The actual dynamics of the observables emerge after 
ensemble averaging. As long as the survival probability 
is not too close to 1/2, a fair representation of P 0 o(i) is 
obtained with about N = 100 realizations as shown in 
Fig. [2] Strongly damped cases evidence the typical fluc¬ 
tuations of random numbers where the observables have 
a precision of 1/y/N. In these cases, individual systems 
maintain a substantial oscillation whose slopes can be 
strongly discontinuous. Thus to obtain coincidence with 
the exact theoretical values within the graphical resolu¬ 
tion (say 1%) one needs about N = 10 000 realizations. 
We tested a binary and Gaussian phase drift distributions 
and the fluctuations and their influence on the precision 
of the ensemble averages persist. However, this is hardly 
a limitation if one is still far from the asymptotic values 
or when one addresses global observables. 


VI. QUANTUM JUMPS VS. QUANTUM DRIFT 
IN A MANY-SPIN DYNAMICS 

Here, we will further assess the differences between the 
QJ and QD in a situation of actual experimental rele¬ 
vance where the density matrix approach is clearly re¬ 
strictive: that of many-spin dynamics. We will address 
the decoherent dynamics of this problem, which is non¬ 
trivial in terms of the density matrix. We are interested 
in the dynamics of a local spin excitation in a system of 
M interacting spins 1/2. Let us say that the state at 
t = 0 is given by the density matrix: 

p 0 = 1 L(i + 2S!), ( 21 ) 

which describes that spin 1 is up-polarized. At very high 
temperature the other spins are not polarized at all, i.e. 
tr[Sfpo] = \ and tr[S?po\ = 0 Vz 7 ^ 1. In order to be 
more specific, let us consider the particular case of a lin¬ 
ear chain with M = 5 that was addressed theoretically 
and experimentally in NMR and where a decoherent cal¬ 
culation is lacking R~il - fT3i| . There, the effective Hamil¬ 
tonian is reduced to nearest-neighbors planar (or XY) 
interactions. The Hamiltonian, using the spin lowering 


S i and raising S/~ operators is: 

M—l 

Hchain = ^ S i+1 + S^ +1 S i ) 

i=1 
M 

+ Y / hn l (s+s~ ( 22 ) 

»=1 

The first term is the XY Hamiltonian, accounting for the 
couplings to nearest-neighbours. We will take Ji.i+i = J 
as the unit of energy. The second one is the Zee- 
man Hamiltonian, where the precession frequencies are 
fli = Wo, and w 0 is the Larmor frequency in the external 
magnetic field. As predicted by the coherent calculation, 
the local excitation />o was seen to propagate as a spin 
wave through the molecular chain and returns to the ini¬ 
tial site in the form of a mesoscopic echo (ME). It is pre¬ 
cisely the wave packet behavior that makes these systems 
promising as quantum channels [ 3 L [~lil - |~i~6| . However, the 
experiments show that these spin waves decohere and at¬ 
tenuate as time pass by. Thus, here we consider that each 
spin is perturbed by a local environment that acts as a 
fluctuating Zeeman field. Thus Qj = uio + 6u>i, where SiOi 
fluctuates with time. 

We will solve this problem by resorting the Jordan- 
Wigner transformation [47] which in this case looks 
5+ 4—> c+ S~ 4—■> c t and S? 4—> ctg. - 1/2 . Thus, 
this many-spin system can be reduced to a one-body 
problem. This transformation remains valid when con¬ 
sidering the random fluctuations in the Zeeman field. 
Then, it is clear that these produce local energy fluc¬ 
tuations Hduji that show up as decoherence in the spin 
dynamics. Note that the fluctuations are naturally de¬ 
scribed within the QD prescription and, conversely, the 
random phases of QD have a direct physical meaning in 
this problem. 

In Figs. [3Ta) and E2A) we compare the coherent and 
decoherent evolutions of the spin wave. The local excita¬ 
tion travels from site 1 to the edge of the chain and there 
it is reflected and returns to the initial spin as a ME. 
We use a decoherence time of about 3 h/J, consistent 
with the experimental observation. j4lj In Figs. [3][c) and 
[3]7I) we show the local density for a single realization for 
both the QD and QJ methods. The QD has a smoother 
profile whereas the QJ resembles the coherent evolution 
until the jump is produced. 

In similar problems, the master equation for the den¬ 
sity matrix could be used to obtain the decoherent dy¬ 
namics [48j . However, with regard to numerical calcula¬ 
tions, it is very demanding. Thus, the QJ method, which 
is based on a stochastic wave function, has proved to be 
more convenient in addressing large systems. [l6j Indeed, 
it is faster and less demanding to perform an average over 
many realizations of the QJ than computing the whole 
density matrix. Since our QD method is also based in 
a stochastic wave function, for large enough systems the 
QD must become more convenient than a density matrix 
approach. Thus, we shall compare the QD with the QJ. 
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FIG. 3. (Color online) Density as function of time at each 
spin for (a) the coherent evolution, ( b ) an average evolution 
with a finite coherence time T co h = 3h/J, (c) a single realiza¬ 
tion of the QD model, and ( d ) a single realization of the QJ 
model. When the spin wave reaches the edge of the chain, it 
is reflected constituting the mesoscopic echo. Note that the 
curves for QD are smooth in density whereas those for QJ 
present a sudden change when the jump is produced. 



FIG. 4. (Color online) The standard error a times vN s as 
function of the evolution time, for N s = 10,100, and 500. 
1 /2 

Note that aN s does not depends on N s but in the QD or 
QJ method. Note that, for the same N s , cr® J /a® D ~ 2. 


Here, T is the evolution time and pu(t) = |z) (*| is the 
mean density at time t at site i obtained from an average 
over N s = 10 5 realizations, where both QJ and QD con¬ 
verge to the the same average dynamics with negligible 
error. In Fig. [I] we show aN for different values of 
N s as function of the evolution time T. Thus, this re¬ 
sult evidences the general tendency of the QD method to 
present a smaller deviation from the mean than the QJ. 
Indeed, for each N s . <jQ J is greater than a^ D by a factor 
of almost 2. Thus, to converge with a given cr, QD needs 
to perform roughly the half of realizations than QJ. 


VII. LOSCHMIDT ECHO 


We will look at the convergence to the average deco- 
herent dynamics for both QJ and QD methods. Due to 
the discontinuous nature of QJ, we expect that the dif¬ 
ference with the mean dynamics is greater for QJ than 
for QD (see Figs. [31(c) and EKd)). This can be quantified 
in terms of the time average of the standard error, 


a (QD,QJ) 






(t) Pi,i 


2 


(23) 


Following the logic of the two previous sections, one 
would be tempted to assign the meaning of decoherence 
to the decay of the oscillations. However, it has been 
clarified that a way to filter away the intrinsic dynamics 
from the observable is by using the Loschmidt echo (LE) 
[33}. This is the amount of excitation recovered after a 
time-reversal procedure implemented in presence of an 
environment. The advantage is that the LE codifies in a 
local observable the losses of non-local correlations. As 
in NMR experiments, this consists in changing the sign 
of the acting Hamiltonian Hq —> —Hq. By using the 
Trotter-Suzuki expansion, the LE can be defined as: 
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where N t = i/di, and E„ is the perturbation’s Hamil¬ 
tonian. Note that, whereas the sign of the Hamiltonian 
changes, the perturbation remains with the same sign. 

In Fig. Ma), we show the survival probability Poo(t) 
and the Losclimidt echo Moo(i) in the underdamped 
regime as function of the time of interaction with the en¬ 
vironment, t. Surprisingly, in the underdamped regime, 
Moo is not a simple exponential but has plateaus when¬ 
ever the reversal starts while the system is at | 0 ) or 
|1). On the contrary, Moo suffers a maximal decay if 
the reversal starts when the system is at a superposi¬ 
tion state [|0> ± 11)] /\/2. When the density is placed on 
one site, decoherent interactions act as a change in the 
global phase, which does not destroy the phase correla¬ 
tions between |0) and |1). Notice that the homogeneous 
exponential decay of the survival probability does not dis¬ 
criminate on the initial state, while the LE does. In any 
case, if one should define a overall rate 'Jc/jle = I/tqle 
from the LE, it would coincide with that observed in 
the oscillation decay, i.e., j^le — = 1/2(tse)~ 1 ■ 

Thus, the LE gives a rationale to the 1/2 factor: decoher¬ 
ent processes are effective on one half of the dynamical 
cycle. 

The overdamped regime is shown in Fig. 0(6) and in 
a log scale in Fig. [5jc) . This last clarifies the differ¬ 
ent decay rates and the difficulties to obtain probabil¬ 
ities around 1/2 from a limited number of realizations. 
The LE has a wider plateau than the survival probabil¬ 
ity at short times. For the same period of action of the 
decoherent processes, the LE signal is higher than the 
survival probability. This fact is consistent with that, 
in order for decoherent processes to be effective, the dy¬ 
namics must first build up the superposition state. The 
decay rate of the LE does not fit with Eq. [20] By using a 
single exponential, the decay rate is I/t^le = 1 . 12 V/h, 
which is slightly smaller than the minimum rate of Pooj 
72 = — ( hjj 0 ) 2 = 1.46V. This indicates that 

the LE gives more weight to the less correlated short-time 
processes where the strong interaction with the environ¬ 
ment does not allow us to create the correlations and 
thus, it should have a slower decay. 


VIII. FINAL DISCUSSION 


We developed and implemented a stochastic model to 
include decoherent processes in quantum dynamics. In¬ 
spired in the Biittiker-D’Amato-Pastawski description of 



FIG. 5. (Color online) Survival probability (thick blue line) 
and the Loschmidt echo (thin black line) as function of the 
time of action of the environment for: (a) r]/V = 0.1 and ( b ) 
T)/V = 2 . 1 , averaged over Ns = 10000 realizations, (c): the 
same as in ( b ) but in log scale. The rate of decay of the LE, 
1 /tse = 1.12V/h, is smaller than the minimum decay rate of 
the Poo, l/i ~2 = 1A6V/H. 


decoherent transport, it boils down to a wave function 
that undergoes smooth stochastic drifts in a local basis. 

Unlike the quantum jumps (QJ) approach, no collapses 
of the wave function occur and phase shifts are intro¬ 
duced in a unitary dynamics. Thus, in QD jumps can 
only appear in the correlations functions, not in the lo¬ 
cal densities. Being an appealing conceptual framework, 
with a clear physical meaning, our QD model results par¬ 
ticularly adapted to deal with extended system. Besides, 
it admits further extensions, ranging from the evaluation 
of currents in transport set-ups to the representation of 
specific many-body interactions. 

Using numerical calculations, we proved that our dy¬ 
namical model is in a full agreement with the decolierent- 
steady-state conductances through the resonant state | 0 ) 
of a decoherent quantum dot, even resorting to a quite re¬ 
stricted ensemble average. For steady-state transport in 
extended systems, a QD evaluation of the wave function 
is, by construction, more efficient than density matrix 
approaches [48| . 

In spin systems, the physical foundation of the QD 
model becomes evident. Decoherence associated to the 
fluctuation of the local energy is a natural ingredient as¬ 
sociated with the fast fluctuation of the local Zeeman 
fields. Thus, we tested the dynamical properties of QD 
model by applying it to a two-spin system and a five 
spin system in presence of decoherence. The first is a 
two-level system (TLS) that oscillates among |0) = |ti) 
and |1) = |4_t). There, a non-trivial quantum dynam¬ 
ical phase transition shows up, which was observed in 
NMR and obtained from the solution of the generalized 
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Landauer-Biittiker equations [29|, [3Cj. We recovered not 
only the exponential damping of the oscillations at low 
rates of interaction with the environment but also the bi¬ 
furcation of the decoherence rates at a critical interaction 
strength. The evaluation of the decolrerent dynamics of 
a five-spin system is also done in connection with NMR 
experiments [4l| . By using rates consistent with the ex¬ 
periment, we show the robustness of the mesoscopic echo 
under decoherence. We also show that, for a given tol¬ 
erance error for an observable, the QD method demands 
about half the realizations than QJ. 

By evaluating decoherence in the TLS through 
Loschmidt echo (LE), we found that the pure states |0) 
= |||) and |1) = ||t) are quite robust against local per¬ 
turbations of the environment. In contrast, the LE, and 
hence coherence, decays faster when the system is in a 
superposition state (|||) ± lit)) /\/2- These results are 
in agreement with the generalJjrend recently observed in 
spin systems through NMR (49| . 


In summary, we proposed a QD model that provides 
a stochastic unitary dynamics of the wave function. Ob¬ 
servable evaluation of observables through QD and QJ 
are naturally parallelizable and thus they result in be¬ 
ing more scalable than density matrix methods [l6j]. This 
quality, which adds to the intrinsic physical significance, 
should make the QD method a suitable tool to address 
dynamical observables in both extended one-body and 
many-body systems. 
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